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Abstract. Symbiotic stars are eccellent laboratories to study a broad range of poorly 
understood physical processes, such as mass loss of red giants, accretion onto compact 
objects, and evolution of nova-like outbursts. As their evolution is strongly influenced 
by the mass transfer episodes, understanding the history of these systems requires fore- 
most to determine which process is at play: Roche lobe overflow, stellar wind accretion, 
or some more complex mixture of both. We report here an interferometric study of the 
symbiotic system SS Leporis, performed with the unique PIONIER instrument. By de- 
termining the binary orbit and revisiting the parameters of the two stars, we show that 
the giant does not fill its Roche lobe, and that the mass transfer most likely occurs via 
the accretion of an important part of the giant's wind. 



1. Introduction 



SS Leporis (HD 41511) contains an evolved M6 III giant and an oversized Al V star 
(see Tab.[T]for a list of its known characteristics), and presents the so-called Algol para- 
dox, as the most evolved star is the least massive of the two, indicating mass reversal, 
through mass loss and mass accretion. The possibility that the system is indeed under- 
going mass transfer is backed up by observations of regular outburst s (|Struve|1930 ) as 
well as from the UV activity from the shell surroundin g the A star (|Polidan & Shore 
1993| ), while the presence of a circumbinary dusty disc ( |Jura et al.|2001 ) hints at the fact 
that the process is not conservative. In this sense, the recent observations of |VerhoeIst 
et al.| ( |200"7] ), indicating that the M star potentially fills its Roche lobe, is of great value. 
We report here new interferometric observations in H and K, obtained with the 



AMBER and PIONIER ( [Berger et aL]|2010| ) instruments attached to the Very Large 
Telescope Interferometer (VLTI). These observations allowed us to directly probe the 
most inner parts of SS Lep, and to unambiguously determine its morphology at different 
epochs, giving rise to a refreshing view of this system. In this contribution, we focus 
our study on the binary system, determining its orbit as well as the M star diameter. We 
will finally show that the giant does not fill its Roche lobe, in a strict sense. 



2. Observations, data reduction and modeling 

Data were collected at the VLTI with the three-telescope spectro-interferometer AM- 
BER and the four-telescope visitor instrument PIONIER. Note that thanks to its addi- 
tional telescope, PIONIER samples the spatial frequencies twice as fast as AMBER. 



1 



2 



N. Blind et al. 





System 




Previous 




This work 


d [pc] 


279 ± 24 


[1] 








P [days] 


260.3 ± 1.8 


[2] 








f(m) 


0.261 +0.005 


[2] 








e 






0.024 ±0.005 


[2] 


0.005 +0.003 


i 






30° ± 10° 


[3] 


143.7° +0.5° 


6a [mas] 






0.53 ± 0.02 


[3] 




6m [mas] 






3.11+0.32 


[3] 


2.296 + 0.007 


M A [Mq] 






2-3 


[2] 


2.71 ±0.27 


M M [M©] 






0.4- 1 


[2] 


1.30 ±0.33 


l/q = M A /M„ 






4 + 1 


[2] 


2.17 ±0.35 



Table 1. Relevant parameters of the SS Lep system, as previously determined 
or from this work: d is the distance, P the orbital period, e the eccentricity, i the 
inclination, and f(m) the mass f unction. For the star s, 6 is the apparent diameter, 
and M the mass. Ref erences: [l] |Van Leeuwen|2007[ [2] |Welty & Wade| 19951 [3] 
IVerhoelst et al.120071 



AMBER data were obtained during 4 different nights in a period of 200 days (more 
than half an orbital period) between November 2008 and April 2009. PIONIER data 
were obtained between October and December 2010 during the commissioning runs of 
the instrument. AMBER observations cover the H- and AT-bands, while PIONIER only 
covers the if -band. For both, we have access to a low spectral resolution R ~ 40 which 
multiplies the number of spatial frequencies sampled. This brings a wealth of infor- 
mation for the parametric modeling, especially for AMBER for which the (w, v)-plane 
coverages were relatively poor (i.e. the number of spatial frequencies sampled is low). 
Typical (w, v)-plane coverages for AMBER and PIONIER can be seen in Fig.[T} 

The data clearly show that SS Lep is a spatially resolved binary whose M giant 
is also resolved in all observations and can be modeled as a uniform disc, whereas the 
circumbinary material is modeled as a gaussian envelope. We tried to detect a possible 
tidal distortion of the giant or matter escaping from its atmosphere by modeling it with 




baseline (m) 



ine (m) 



Figure 1 . Typical (u, v)-plane coverage of SS Lep as obtained with AMBER (left) 
and PIONIER (right). Each point represents a spatial frequency sampled by the 
instrument. The concentric dotted circles represent projected baselines lengths equal 
to 15, 33, 50, 65, 80 and 1 15 metres 
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Figure 2. Model independent image reconstructions of SS Lep obtained during 
the four runs with PIONIER. The resolved M giant (about 2 mas large) and the A 
star are clearly identified, and separated by roughly 5 mas. The distortion of the giant 
in the images is most certainly due to an asymmetric PSF rather than to a real tidal 
effect. A few faint artefacts are visible on the periphery of the images. 



an elongated uniform disc. Results were not conclusive as we lack the longest baselines 
required to measure distortions of the order of a few percents. The spatial resolution was 
not sufficient to resolve the putative shell or an accretion disc around the A star. Given 
the longest baselines of our observations (130 m with a spatial resolution around 1 mas), 
this agrees with the 0.5 mas size estimated from the Spectral Energy Distribution (SED) 
in |Verhoelst et al.| ( |2007| ). The model we used to fit the interferometric data therefore 
comprises six degrees of freedom: the relative flux contribution of two components 
of the system; the binary separation and its orientation; the size of the M giant and the 
size of the circumbinary envelope. The lack of short baseline information prevents us to 
constrain the spatially extended emission due to the large circumbinary disc (FWHM ~ 
12 mas). We therefore focus on the binary in the following. 



3. Image synthesis 

With its four telescopes, PIONIER allowed us to perform a reliable model independent 
image reconstruction of SS Lep with the MIRA software (Thie baut||2008| ) for the four 
observations (see Fig. [2]). Each image clearly shows the binary nature of SS Lep, the 
separation being slightly smaller than 5 mas. From one observation to the other, we can 
observe the rotation of the system. 



The A star and its shell have an expected spatial extension of 0.5 mas ( Verhoelst 



et al.||2007| ) so that we do not expect to resolve them with the VLTI baselines. There 



fore, the size of the spot corresponding to the A star defines more or less the point 
spread function (PSF) of the image, about 1 mas. The M giant being the most lumi- 
nous component of the system in the //-band, we identify it on images as the darkest 
spot. With respect to the A star, we clearly see that it is spatially resolved and measures 
approximately 2 mas in diameter, while the orbit size is roughly 5 mas. We expect the 
distortion observed on the images to be due to an asymmetric filling of the (w, v)-plane 
(implying a non-circular PSF on the reconstructed image) rather than to a real tidal 
distortion. As a matter of fact, the tidal distortion would be around 7%, i.e. less than 
seen in the images. Additionally, its orientation in the images corresponds well with 
the asymmetry observed in the corresponding (w, v)-planes. It was actually not possible 
to image the circumbinary disc because of the lack of data with short baselines. 
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Figure 3. SS Lep's best orbit (dashed line) obtained by combining previous radial 
velocities ( Welty & Wade|1 995 ) with our astrometric measurements. The A star is 
indicated as the dot in the centre. AMBER and PIONIER points are represented by 
the triangles and crosses, respectively, with the error bars. The corresponding points 
on the best orbit are the closest circles from each point. 



4. The binary 

4.1. The orbit of SS Lep 

To compute the most reliable orbit possible, we combined the eight astrometric posi- 
tions of the binary (separation and orientation of the stars) obtained from our interfer- 
ometric observations, with the radial velocities of [Welty & Wade| ( [l995| ). The best-fit 
orbit is shown in Fig. [3j and the derived orbital parameters are listed in Tab. [2] The 
errors on the orbital elements are estimated via Monte-Carlo simulations. The incli- 
nation angle of 143.7° and the close to zero eccentricity are in agreement, but more 
precise, than the previous estimates by |Welty & Wade| This result definitely invali- 
dates the periastron-passage mass transfer scenario of |Cowley| ( |1967| ), which required 
an important eccentricity (e = 0.134) to explain the regular outbursts of the system. 



Semi major 
axis a 


Inclination 
angle i 


Eccentricity 

e 


Longitude of the 
ascending node Q 


Argument of 
periapsis co 


4.492 + 0.014 mas 


143.7 ± 0.5° 


0.005 ± 0.003 


162.2 ± 0.7° 


118 + 30° 



Table 2. Orbital parameters of SS Lep obtained by combining previous radial 
velocities ( Welty & Wade|1995 ) with our 8 astrometric measurements. 



4.2. The distance and the masses 

Combining the orbital parameters with the binary mass function estimated by | Welty & 
|Wade| ( [T9^5l ), we can estimate the individual mass of the stars, and, thus, the mass ratio. 
The main source of uncertainty in this estimate resides in the distance, as determined 
by Hipparcos. Using the distance and the angular separation of the two stars, we obtain 
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Figure 4. Representation of the Roche equipotential (solid line) for a mass ratio 
1 lq = 2.2. The limb-darkened diameter of the M giant is in dashed line, and the A 
star one is the dark dot. 



a - 1.26 ± 0.06 AU, and thus, through Kepler's third law, the total mass of the system 
is estimated as 4.01 ± 0.60 M Q . Thanks to the inclination that we determined and to the 
binary mass function, we can derive the individual masses and the mass ratio: Ma = 
2.71 ± 0.27 Mo, M M = 1.30 ± 0.33 M Q , and l/q = M A /M M = 2.17 ± 0.35. 

4.3. The M star diameter 

We measure a uniform disc diameter for the M star 6m,ud = 2.208 + 0.012 mas. The 



previous VINCI observations of Verhoelst et al. led to a higher value of 2.94 + 0.3 mas, 



most likely because theirs was the result of a one-year survey of the source, without 
any phase information in the interferometric data. This involved to model the system 
as a symmetric object, and to mix up the interferometric signatures of the rotating 
binary with the resolved giant one. Depending on the authors ( |Hanbury Brown et~aT 



T9741 |Davis et al.] [20001 °r |Verhoelst et~^|2007] ), the conversion factor to a limb- 



darkened diameter differs by a few percents and we fix it to 1.04. This leads to a 
limb-darkened diameter equal to 6m,ld = 2.296 + 0.013 mas. Our results is still in very 
good agreement with the limb-darkened diameter estimated from the SED of |Verhoe"Ist 



et al. (6m,ld = 2.66 + 0.33 mas). Taking into account the uncertainty on the distance, 
the M giant radius is thus Rm = 69 ± 3 R Q (almost twice as small as previously thought), 
and the surface gravity is therefore logg ~ 0.8. 



5. The mass transfer process 

The circumbinary disc and the envelope around the A star, as well as the activity seen in 
the ultraviolet, are good hints for the occurrence of mass transfer. The disc additionally 
indicates that the process is not conservative. With the parameters derived above, we 
are able to show that, conversely to previous results, the M giant only fills about 85+3% 
of its Roche lobe (Fig. [4]), and the mass transfer is more likely due to the accretion of 
the M giant wind. If we assume the star's surface to follow an equipotential surface, 
we expect a tidal distortion of the giant from 5 to 7%. The wind speed being lower 
than the orbital one (v or b = 48 kms" 1 ) the system is possibly in the particular case of a 
wind Roche lobe overflow ( |Mohamed & Podsiadlowski|2007| ), where a substantial part 



of the stellar wind can be accreted. From the simulations of |Nagae"eT al. ( 2004]) the 



wind accretion efficiency for these velocities is around 10% for a mass ratio equal to 
1. The accretor in SS Lep being the most massive component, the efficiency could be 
significantly higher. However we do not expect that all the wind can be accreted because 
the circumbinary disc is a proof of the non conservative behaviour of the mass transfer 
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process. Assuming the M giant looses 10" 7 M yr _1 ( Dupree||l986| ) through wind, the 
accretion rate onto the A star should be between 10" 8 and 10" 1 M yr _1 . It is thinkable, 
however, that some enhanced mass loss takes place in such kind of systems, and this 
is possibly only a lower limit. This is still much lower than the value of 10" 4 M yr _1 
quoted by |Verhoelst et al.| ( |2007] ). 



6. Conclusion and future work 

We have presented here the results of our observations, and we focused on the binary. 
After having computed the characteristics of the orbit, we demonstrated that the mass 
ratio is lower than previously thought and that the M giant does not fill its Roche lobe. 
However the system is in a configuration where a substantial part of the giant's stellar 
wind can be accreted by the A star. The refreshing and precise view of SS Leporis 
given by these observations paves the way for the study of interacting binaries with 
interferometry, in combination with spectroscopic studies. 

The interferometric data also gave us information about the relative luminosity of 
the three components of the system (the M and A stars, and the circumbinary disc). 
There is still some work to fully exploit these data. In particular, we lack a low resolu- 
tion spectrum of SS Lep between 1.6 and 2.5 fim to compute the absolute luminosity of 
each component and extract more specific information. The current data also present 
good hints of unmodeled material escaping the system. We plan to observe SS Lep with 
the full high spatial resolution capacity of the VLT, with NACO and PIONIER to study 
in more details the binary and the circumbinary disc, and extract as much information 
as possible about the mass transfer and the interactions. We also plan to use AMBER 
to observe the M giant wind and detect an asymmetry due to the high efficiency of its 
accretion onto the A star. 
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